Absolute determination of Ds branching ratios 
and fi^^ extraction at a neutrino factory 



G. De Lellis^), P. Migliozzi^) and P. Zucchelli^) 

1) Universita Federico II and INFN, Napoli, Italy 
2) INFN, Napoli, Italy 
3) CERN, Geneva, Switzerland and INFN, Ferrara, Italy 

Abstract 

A method for a direct measurement of the exclusive Ds branching ra- 
tios and of the decay constant fo, with a systematical error better than 
5% is presented. The approach is based on the peculiar vertex topol- 
ogy of the anti-neutrino induced diffractive charm events. The statistical 
accuracy achievable with a neutrino factory is estimated. 



1 Introduction 



The experimental knowledge on leptonic decays is rather poor. Currently, 
the branching ratios for Dg — > Ih' decays are estimated by the Particle Data 
Group II to be BR{Ds -> fJ-iy) = (4.6 ± 1.9) x lO'^ and BR{D, tv) = 
(7 ± 4) X 10^^. These large uncertainties translate into a large uncertainty on 
the extraction of the decay constant f^^ . 

In this paper we propose a method to build an almost pure sample of 
from diffractive events which allows the extraction of most of the Ds branch- 
ing ratios and in particular of purely leptonic decays from which fjj^ can be 
extracted with a systematic error better than 5%. A statistics capable of ex- 
ploiting this systematic error can be accumulated at neutrino factories. Once 
fiy^ is measured with such an accuracy, one would feel more confident about 
extrapolating to the decay constants in the B system, fs and /b, , which are 
crucial quantities for a quantitative understanding of B^^j ~ (s) oscillations and 
the possible extraction of the CKM matrix elements Vtd or Vts 

This paper is organised as follows: in Section ^ we discuss the available data 
on neutrino and anti- neutrino induced diffractive D^*' [^production. In Section || 
the leptonic decay mechanism is discussed and the available experimental 
determinations of fn^ are reviewed. A method for a direct evaluation of Dg 
branching ratios and fo^ measurement is discussed in Section ^ together with 
the evaluation of the accuracy achievable at a neutrino factory. In the last 
section we give our conclusions. 



2 (Anti-)Neutrino induced diffractive Ds pro- 
duction 

In charged-current (CC) interactions, the W boson can fluctuate into a charmed 
meson. The on-shell meson is produced by scattering off a nucleon or a nucleus 
without breaking up the recoiling partner. The diffractive production mecha- 
nism (Ufj^N ^~D*~^N) is shown schematically in Fig. ^ The same mechanism 

applies to D^*^ production, but the oi*^ one is Cabibbo favoured by a factor 



This process has been observed in previous experiments ^, ^. In par- 
ticular, one of these experiments Q has shown the evidence of a Ds diffractive 
production through the direct observation in nuclear emulsion of the decay chain 
D*+ ^ D+-f, D+ ^ T+Ur. T+ ^ ^^+yry^.■ 

Unlike deep-inelastic cross-sections, the diffractive cross-sections are pre- 
dicted to be the same for v and v and for proton and neutron targets. The 
observed neutrino and anti-neutrino induced diffractive Ds rates relative to 
CC interactions on an isoscalar target have been measured to be (1.5 ± 0.5) x 
10~^ and (2.6 ± 0.9) x 10"'', respectively. This is compatible with a 1 : 2 

^In the following D^*^ stands for either Ds or D*. The same notation is also used for D^*\ 
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ratio, as implied by the equality of difFractive cross-sections and the ratio of 
inclusive cross-sections S. It is worth stressing that these numbers have been 
obtained searching for particular Dg decay channels in Ds diffractive produc- 
tion: therefore the corresponding branching ratio has to be known to get the 
absolute production rates Q. 

The energy dependence of the difFractive cross-section has also been investi- 
gated. In the neutrino energy intervals 10 30, 30 50 and 50 ~ 200 GeV the 
observed rate per CC is (1.8±0.7) x 10-^ (1.3±0.6) x 10^^ and (1.6±0.7) x 10"^ 
respectively, so that no variation is detected at the available statistical level Q . 
Theoretically, no steep variation of this relative rate in the 10^200 GeV neutrino 
energy interval is expected. 

Taking into account the corresponding branching ratio, the neutrino difFrac- 
tive ni*^ production rate has been evaluated to be (2.8 ± 1.1) x lO^'^/CC 
An independent evaluation of the same production rate has given (3.2 ± 0.6) x 
lO^'^/CC 0. The weighted average of the two evaluations gives a neutrino in- 
duced difFractive production rate oF (3.1 ± 0.5) x 10^"^/CC. Since the neutrino 
to anti-neutrino CC production rate is 2 : 1, the average anti-neutrino induced 
-Oi*-* diffractive production rate is (6.2 ± 1.0) x lO^'^/CC. Therefore, the com- 
bined analysis gives an accuracy of about 15% for the diffractive production 
rate. 

For a detailed theoretical review of lepton induced diffractive production we 
refer to Refs. i | | 0. 




Figure 1: Diagram for neutrino induced diffractive D*^ meson production. 



3 Leptonic decays 

From a theoretical point of view, purely leptonic decays of charged mesons 
are the simplest ones to describe. The effect of the strong interaction can be 
parameterised in terms of just one factor, called the decay constant. Unlike 
semi-leptonic decays, where (and hence the form factors) varies event by 
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event, leptonic decays have a fixed value (q^ = Af^, where Af is the mass of 
the initial meson). 




By neglecting radiative corrections, the decay rate of a charged pseudo-scalar 
meson Mqq to Ivi is 

/ 2 \ ^ 

T{Mq, ^ Ivi) = -^x I Fq, p x/l, X M xmfx[l-^J 

where /m is the decay constant, Vqq is the CKM matrix element, and mi 
and M are the masses of the lepton and charged meson Afgg, respectively. The 
decay constant /a/ is a measurement of the probability amplitude for the quarks 
to have zero separation, which is necessary for them to annihilate. 

By using the previous notation, the Dg leptonic branching ratio can be 
written as 



BR{Ds ^ l-vi) = ^x\VQq\^ xfl^ X TD^ X Md^ x mj x [l - j^j 
where tjj^ is the Dg life-time. 

Decay constants for pseudo-scalar mesons containing a heavy quark have 
been predicted with lattice QCD, QCD sum rules and potential models, but 
due to the non-perturbative character of the calculations they vary significantly 
in the predictions. The predicted value for fn^ lies in the 190 360 MeV 
range [pl]| . 

A method for extracting fo^ is to measure the leptonic decay modes of 
Ds- Due to helicity suppression, only the muonic and tauonic decay modes 
have an appreciable branching ratio. Unfortunately, the measurements of the 
leptonic Dg decays are rather scarce. A summary of the available measurements 
together with the determinations of fjj^ is shown in Table |l|. The errors on fu^ 
are larger than 30% and, despite the scanty statistics available, the systematic 
error dominates. This is mainly due to the uncertainty on the normalisation 
used for the determination of the leptonic decay branching ratios. 

Once the absolute value of is measured with an accuracy of 10% or bet- 
ter, one would feel more confident about the predictions of the decay constants 
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Experiment 


Channel 


/d, (MeV) 


WA75 


12^ 


Ds ~* (I 


232 ± 45 ± 52 


CLEO I Ml 


Ds fi 


344 ± 37 ± 52 ± 42 


CLEO IlH 




280 ± 19 ± 28 ± 34 


E653[ 


15| 


-Ds — * M 


194 ± 35 ± 20 ± 14 


BEATRICE|16| 


Ds —* fi 


323 ± 44 ± 12 ± 34 


BESJ17I 


Ds I 


430+J30 ±40 




Ds^T 


309 ± 58 ± 33 ± 38 



Table 1: Summary of the available experimental determinations of the decay 
constant fo^ ■ 

in the B system, fs and fs^ i which are crucial quantities for a quantitative 
understanding of B'^^^ — B'^^^ oscillations and the extraction of VtdiVts) from 
them. 

In the following we describe a method for the extraction of the Dg leptonic 
branching ratios which minimises the systematical error due to the normalisation 
of the sample. It is based on the direct observation of Ds ' produced in anti- 
neutrino induced diffractive interactions and a topological selection of the events 
which improves the purity of the sample used for the normalisation. 

4 Direct evaluation of Dg branching ratios and 
fjj^ measurement 

In the following we consider a detector with the capability to exploit vertex and 
decay topologies with micron precision. We assume to use an emulsion target 
(see for instanc e |19(| ), although the method we are proposing would work with 
a CCD target as well. This allows the detection of short lived particles 
with path length larger than 10 fim. In order to measure the charge and the 
momentum of hadrons and muons, electromagnetic and hadronic showers, a 
possible detector is the one described in Ref. ||^ once equipped with such an 
emulsion target. 

4.1 Event yield and flux at a neutrino factory 

For this study we assume a muon beam energy = 50 GeV; length of the 
straight section, L = 100 m; muon beam angular divergence, 0.1 x m^/E^, 
being the muon mass; muon beam transverse size ax = Cy — 1.2 mm. 

The target is made of N emulsion bulks 1 m diameter and 10 cm thick 
each. By locating the detector 1 km far from the neutrino source, we obtain the 
neutrino energy spectrum as shown in Fig. |^. 

The use of nuclear emulsions is limited by the interaction overlapping. A 
density of interactions of about 20 cm"^ is affordable. Therefore, aiming at a 
statistics of 10^ z/^ interactions we assume N — 16. For 1 year data taking, a 



5 



neutrino flux of = 1.6 x \{f /lQ{v/ s) through the detector would be required. 
It corresponds to about 3.5 x 10^'' decaying muons. Under these assumptions we 
obtain an overall density of interactions of about 16 per cm^, which is reasonable. 

4.2 Kinematical selection: deep-inelastic versus diffrac- 
tive events 

Deep-inelastic charmed events are kinematically quite different with respect to 
diffractive ones. In particular, the kinematical variable i, defined in Fig. ^, al- 
lows a very high rejection power against the background with a good efficiency 
of the signal. The t determination relies on the momentum determination of 
the charmed particle. In the E531 emulsion experiment a resolution 

in measuring the charmed hadron momentum better than 15% was achieved 
by using a likelihood technique. In the following we make the conservative 
assumption of momentum resolution in measuring the charmed hadron momen- 
tum of 30%, 50% and 100%. The reconstructed t distributions for diffractive 
and deep-inelastic events are shown in Fig. |^. 

In order to select diffractive events we apply the cut i < 1.7 GeV. Cuts on the 
flight length (> 10 /im) and on the kink angle {Okink > 15 mrad) are also applied. 
The efficiency for the signal and the correspondent background are reported in 
Table || From this table we can see that, independently of the momentum 
resolution, the fraction of deep-inelastic charm events with diffractive topology 
surviving the kinematical cut is Sdis 0.4% , while the signal efSciency ranges 
from about 50 to 80%. In the following we assume a momentum resolution of 
50%. 



Ap/p 


£*s(%) 




30% 


0.3 ±0.1 


83.3±0.8 


50% 


0.4 ±0.2 


72.5±1.0 


100% 


0.4 ±0.2 


51.8±1.1 



Table 2: Efflciency for the signal, Ed^, and for the background, edis, assuming 
three different momentum resolutions. 



4.3 Topology of neutrino induced diffractive charm events 
and background 

In the ni*'' or diffractive production only a muon is produced at the in- 
teraction point (primary vertex), besides the charmed meson. Therefore, these 
events are characterised by a peculiar topology: two charged tracks at the pri- 
mary vertex, one of them being a short lived particle. 

Note that the excited states produced undergo the following decays without 
leaving an observable track: 
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D*-^D-j Si? = (94.2 ± 2.5)% g 

D*-^D-7:" Bi? = (5.8±2.5)% § 

D*- ^ D^ir- Si? = (67.7 ± 0.5)% g 

D*-^D-7:" 5i? = (30.7 ± 0.5)% § 

D*- ^ D--f Si? = (1.6 ±0.4)% 11. 

A source of irreducible background is the diffractive production of D^*^' 
Its contamination relatively to the signal is given by 



£_D(.) 



Vcd ' 



X [riD- + Vd- X BR{D*- 



where BR{D*- D-) ^ 0.323 ± 0.006 §, 77^,- and r]^.- are the fractions of 
diffractively produced and D*^ respectively. We assume that these fractions 
are the same for D~ as for D~ {rijj- = rjjj,- — 0.5). The latter can be extracted 
by the NuTeV results: (j{v^N ^i^DsN) = (1.4±0.3)fb/nucleon and cr(i/^7V 
H-D*N) = (1.6 ± 0.4)fb/nucleon g. Finally we get eo- = (3.3 ± 0.8)%. 

Neutrino induced quasi-elastic charm events are characterised by the topolo- 
gies shown in Fig. |6| (see Ref. |Q). Therefore, they are similar to the diffractive 
ones, but with a cross-section twice as large. Since anti-neutrinos cannot induce 
quasi-elastic charm production while diffractive production is the same for both 
V and V, our interest is focused on 9 induced events. Therefore we make the as- 
sumption that all the events with the above topology are due to Ds diffractive 
production. 

The charm production in D deep-inelastic interactions and the event topol- 
ogy have been studied by using the HERWIG event generator an event 
generator based on JETSET Q and LEPTO Q and the energy dependence 
of the charmed fractions reported in Ref. | |2^p|. The average charmed fractions, 
convoluted with the anti-neutrino spectrum are F^o = 61%, F^- = 26%, 
Fj-^- = 7.3% and F. - — 5.7%. The kinematics signal has been modeled accord- 
ing to Refs. |6[ 1^, g, ^, |lO| by using an event generator developed within the 
CHORUS Collaboration |5 1 . 

The contamination of the diffractive sample from deep-inelastic events can 
be written as 



aiPf^N ijL+CX) 1 



Efake = , X:5-X(^g- + ^A" ) ^ ffake^^SMs (1) 



where 



^We assume that the charmed fractions are equal for both 1/ and u as implemented in the 
event generators. 
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fake/Sfak£ 


£fake{%) 


15% 


22% 


0.037±0.009 


30% 


34% 


0.04±0.01 


50% 


53% 


0.04±0.02 


100% 


101% 


0.04±0.04 



Table 3: The relative and absolute error on Efake as a function of the relative 
error on TZ. 



We take the charm production in v interaction to be 3% . By using the 
charmed fractions given above, (Fd- +F^-) = 31.7%. The factor ffake = (6.0± 
0.1)% denotes the fraction of deep-inelastic charmed events faking a diffractive 



topology. £dis gives the efficiency of kinematical cuts as explained in Section 4.2 
As discussed in Section ^ the experimental determination of H has an accu- 
racy of about 15% which gives Efake = (0.037 ± 0.009)%. The Sfake value and 
its error are reported in the third column of Tabic ^ as a function of the error 
on 1Z. 

From the numbers given above it turns out that the little knowledge we 
have about the diffractive charm production cross-section plays a role only in 
the evaluation of the deep-inelastic contamination, namely a term of the sys- 
tematic error. Even if the ratio TZ had an uncertainty of 100%, at Str the relative 
systematic error on the branching ratios would be < 0.16% (see Table ||), so 
that the e^jc.) contribution is dominant. Therefore, the overall systematic un- 
certainty does not depend at all on the 1Z accuracy {egys = (3.3 ± 0.8)%). 



4.4 Description of the method 

An almost pure sample of D'^ from diffractive events, with a small contam- 
ination of and A^ produced in deep-inelastic events and in diffractive 
£)(*)- production, can be built by using diffractive production from anti- 
neutrinos. The normalisation to determine the Ds absolute branching ratios is 
given by the number of events with a vertex topology consistent with one /li plus 
a short lived particle. No model dependent information is used to define the 
normalisation. 

It is also worth noting that, in particular, the contamination of and A^ 
events does not affect the Dg ^ t channel. Indeed, such events would present 
a unique topology with two subsequent kinks. An event with a double kink has 
been recently observed in CHORUS (see Ref. Q). 



4.5 Measurement accuracy at a neutrino factory 

At present there are no experiments with both an adequate spatial resolution 
to fully exploit the diffractive topology and a sufficient anti-neutrino induced 
CC event statistics. Therefore, the method proposed in this paper could only 
be exploited with the above-mentioned detector exposed at a neutrino factory. 
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Channel 


PDG BR [1| 


This method 


Da ^iv 


(4.6 ± 1.9) X 10"^ 


(±0.55±0.15) X 10"^ 


Ds TU 


(7 ±4)% 


(±0.17±0.23)% 


Ds (t>lv 


(2.0 ±0.5)% 


(±0.08±0.07)% 



Table 4: Statistical and systematic accuracy achievable in the determination of 
the Ds absolute branching ratios, assuming a collected statistics of 10'' i?^ CC 
events. The central values are taken from Rcf. ||l[. 

Let us assume to collect 10^ vCG events into an emulsion target and to 
have a detection efficiency of about 73% for the Ds decays (see Table ||). By 
assuming a vertex location efficiency of about 50%^ and assuming a v diffractive 
production rate of 6.2 x 10~''/CC, we expect to detect a number of Ds equal to 
Nd^ = 10^ X 6.2 X 10-3 X 0.73 x 0.5 ~ 2.3 x 10''. 

The expected accuracy on the determination of the Ds branching ratios is 
shown in Table § for a few channels, together with the current status. To 
compute the expected number of events in each decay channel we have used 
the central values (shown in Table ^ together with their errors) given by the 
Particle Data Group Q|. 

By using the relation given in Section ^ and the measured branching ratios 
given in Table ^, the decay constant Jd, can be extracted. If we collect 10^ 
CC interactions we get fn^ ~ 288 ±4 \stat ±5 \sys MeV where the central value 
is taken from Ref. Q. 

5 Conclusions 

Ds branching ratios are affected by large uncertainties, mainly due to the diffi- 
culty of defining a pure Ds starting sample. These uncertainties translate into 
the determination which is in turn affected by large errors. A method for 
the evaluation of the Ds branching ratios and of the fo^ decay constant with a 
systematic accuracy better than 5% has been presented. 

The idea is to build an almost pure sample of -Ds's by means of the anti- 
neutrino induced diffractive Dl ' production. The vertex topology of these 
events is extremely simple; there is only a muon and a short lived charmed 
hadron, the Ds- This peculiarity makes the contamination of D~ and A" from 
anti-neutrino deep-inelastic interactions negligible. The diffractive D^*^ pro- 
duction yields a contamination of about 3% which is the dominant systematic 
uncertainty. 

In order to make the statistical error at the same level as the systematic one, 
a copious (0(10^")) number of anti-neutrino interactions is needed. 

At present there are no experiments with both an adequate spatial resolution 

^This efficiency accounts for the electronic detector reconstruction and the automatic lo- 
cation of the event vertex inside the emulsions. 
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to fully exploit the difFractive topology and a sufRcient anti-neutrino induced 
CC event statistics. Therefore, the measurement could only be performed at a 
neutrino factory. 
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Figure 3: Neutrino energy spectrum in the assumption of Im diameter detector 
located 1 km far from the neutrino source. 
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Figure 4: The 4-momentum transfer f in a neutrino CC interaction with charm 
production, q' is the charmed hadron 4-momentum. 




Figure 5: t distributions for both diffractivc [Ds) and deep- inelastic {fake) inter- 
actions with a diffractive topology. A resolution of 30% (left plot), 50% (middle 
one) and 100% (right one) on the charmed hadron momentum measurement is 
assumed. 
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a) b) c) 



Figure 6: Topology of the quasi-elastic charm neutrino induced events in the 

case of the reaction a) i>^n /i~A+, b) v^n /i~S+(E*"'') and c) i^^p 

The particles inside the box represent the A+ decay products. 
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